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Estimates of hadron azimuthal anisotropy from multiparton 
interactions in proton-proton collisions at J~s = 14 TeV 
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Abstract. We estimate the amount of collective "elliptic flow" expected at mid-rapidity in proton-proton (p-p) colli- 
sions at the CERN Large Hadron Collider (LHC), assuming that any possible azimuthal anisotropy of the produced 
hadrons with respect to the plane of the reaction follows the same overlap-eccentricity and particle-density scalings 
as found in high-energy heavy ion collisions. Using a Glauber eikonal model, we compute the p-p eccentricities, 
transverse areas and particle multiplicities for various phenomenological parametrisations of the proton spatial density. 
For realistic proton transverse profiles, we find integrated elliptic flow V2 parameters below 3% in p-p collisions at 
V^s = 14 TeV. 

PACS. 13.85.-t- 13.85.Hd-25.75.Ld 



1 Introduction 

With increasing energies, hadronic collisions are characterised 
by a larger and larger number of produced particles issuing 
from the fragmentation of a growing number of partons in- 
volved in the reaction. In the extreme case of high-energy nucleus 
nucleus (A-A) reactions, the multiplicity of released partons is 
so large - e.g. ^(lOOO) at midrapidity in central Au-Au colli- 
sions at RHIC energies ( = 200 GeV) - that they interact 
strongly among each other leading to "hydrodynamical flow" 
behaviours [1]. A distinct consequence of such partonic ex- 
pansion effects is an azimuthal anisotropy of the final particles 
produced, dN/dA§, with respect to the reaction planqj, A(j) = 
(j) — <Prp. At RHIC energies, the large "elliptical" anisotropy 
measured in the data has provided detailed information on the 
degree of collectivity reached in the early stages of nuclear col- 
lisions (see e.g. [2] for a recent comprehensive review). 

The parton density inside hadrons grows very rapidly with 
increasing energies, as more and more gluons share smaller and 
smaller fractional momenta x = p par ton/ Phadron ( see e -g- [3, 4] 
and refs. therein). Thus, at high enough energies the proton it- 
self can be considered as a dense and extended partonic ob- 
ject, and proton-proton collisions can be viewed very much 
like collisions of "light nuclei" composed of their constituent 
gluons. The possibility of having multiple parton interactions 



1 The reaction plane of a collision is spanned by the vector of the 
impact parameter b between the centers of the two colliding objects, 
and the beam direction. Its azimuth is given by $>rp . 



(MPI) happening simultaneously at different impact parame- 
ters in hadronic collisions has been discussed since long in 
the literature [5]. Signals of multiple parton scatterings have 
been observed experimentally in p-p collisions at centre-of- 
mass (cm.) energies i£?(l TeV) in charged particle multiplic- 
ity distributions [6] as well as in multi-jets events [7-9]. At 
LHC energies, recent developments of general-purpose Monte 
Carlo (MC) event generators such as PYTHIA [10] and HER- 
WIG [11] include an impact-parameter-space description of the 
p-p collision to account for MPI effects. As a result of the large 
number of low-x gluons and increasing multi-parton collisions, 
the MC predictions for the total number of produced particles 
in "central" p-p collisions at y/s = 14 TeV, are pretty large - 
up to five hundreds hadrons - similar to those measured e.g. in 
intermediate-size nuclear (Cu-Cu) reactions at -Js^ = 200 GeV 
[12] where a significant elliptic flow has been observed in the 
data [13]. Given this situation, it does not seem unjustified to 
contemplate the possibility of having some type "collective" 
behaviour also in the final-state of proton-proton collisions at 
LHC energies [14]. 



In this paper, we explore the possibility of observing gen- 
uine collective expansion effects, such as an azimuthal anisotropy 
with respect to the plane of the reaction, in p-p collisions at 
yfs = 14 TeV. We note that the partons involved in the scatter- 
ings of relevance for such collective effects have rather moder- 
ate virtualities of the order of the average transverse momen- 
tum expected in a minimum-bias p-p collision at the LHC, 
i.e. <??(0.7 GeV). Recent attempts to estimate similar effects 



2 



D. d'Enterria et ah; Estimates of collective hadron azimuthal anisotropy in p-p collisions at \fs = 14 TeV 



have been considered in the context of percolation [15], colour- 
dipole saturatiorQ [16] and (viscous) hydrodynamics [17] mod- 
els. In our work, we take a more empirical (i.e. less model- 
dependent) approach and we simply assume that the eccentric- 
ity and multiplicity dependence of the elliptic flow experimen- 
tally observed in high-energy heavy ion collisions (Section [2]i 
holds also for proton collisions at the LHC. Using a geometri- 
cal eikonal approach (Section[3]), we determine the eccentricity 
of the transverse overlap region in p-p collisions for different 
proposed proton spatial densities (Section|4]i, and then we esti- 
mate the resulting associated azimuthal anisotropy for various 
impact-parameters (Section[5]l. 



2 Elliptic flow scaling in ultrarelativistic heavy 
ion collisions 

One of the main results at RHIC is the observation [18, 19] of 
a large harmonic modulation of the azimuthal distribution of 
the produced hadrons, dN /dA§, with respect to the reaction- 
plane in A-A collisions with non-zero impact parameter, i.e. 
with a lens-shaped overlap zone between the colliding nuclei. 
The measured preferential "in-plane" emission is consistent with 
an efficient translation of the initial coordinate-space anisotropy 
into a final "elliptical" asymmetry in momentum-space through 
rescatterings between the produced partons in the early stages 
of the collision [20, 21]. The strength of the elliptic flow asym- 
metry in A-A collisions is quantified via the second Fourier co- 
efficient v'2 = (coi(2A(j))) of the azimuthal decomposition of 
the single inclusive hadron spectrum with respect to the reac- 
tion plane [22,23] 

d 3 N 1 d 2 N / " \ 
d 3 p 2np T dp T dyy ~ { J 

The large integrate^ V2 ~ 0.06 measured in Au-Au collisions at 
y/s^ = 200 GeV [18, 19] is well reproduced by hydrodynami- 
cal simulations [20,21] that predict the development of collec- 
tive motion along the pressure gradient of the system, which 
is larger for the in-plane directions parallel to the smallest di- 
mension of the lens-shaped overlap zone. Within this scenario, 
one naturally expects V2 to be directly proportional to (i) the 
density of produced particles in the transverse plane, as well as 
to (ii) the original spatial anisotropy of the system quantified 
by its eccentricity, usually defined as e = (j/ 2 — tt) / (x 2 + i/ 2 ) 
(where x,y are the transverse dimensions of the "lens", and 
the average is taken over the initial nuclear profile with some 
weight, see later) [24,25]. 

Experimental heavy ion data at different cm. energies with 
different colliding systems and varying centralities [26] indeed 



indicate a proportionality of the vijz ratio with the particle 
multiplicity normalised by the size of the system, i.e. V2/e °c 
(dN /dy) /A±, where dN /dy is the total hadron multiplicity^ N 
per unit rapidity y, and the overlap transverse area Ax is cal- 
culated via a Glauber geometrical model (see e.g. [27]) with 
standard Woods-Saxon distributions for the initial nuclear dis- 
tributions [28]. An attempt to parameterise the observed V2/E 
dependence on particle density based on an "incomplete equi- 
libration" model has been developed in Refs. [29, 30] in terms 
of the Knudsen number K = X/R, where A, is the mean free 
path of the interacting partons and R a typical length scale (e.g. 
the radius) of the system. By definition, K 1 is the mean num- 
ber of collisions per particle (i.e. the medium opacity), and 
the ideal hydrodynamics limit corresponds to K — > 0. Since 
X = 1/ (a gg p), where a gg is the effective parton-parton (mostly 
gluon-gluon) cross section, and since p(x) = 1/ (iA±) ■ dN /dy 
is the (time-dependent) density of the medium, for a typical 
time x = R/c s (where c s is the medium speed of sound) one 
can write the (inverse of the) Knudsen number as 



2 Strictly speaking, the origin of elliptic flow in this work is of a dif- 
ferent nature (initial parton distributions) than the final-state collective 
effects of interest here. 

3 The V2 coefficient is a function of rapidity and transverse momen- 
tum, and as such it is often referred to as differential flow. Hereafter, 
we will be interested in integrated flow, namely the value of the V2 co- 
efficient averaged over transverse momentum and rapidity in a given 
event. 



G gg c s dN 
Ax ~dy' 



(2) 



Within this model, the dependence of v'2/e on dN /dy in A-A 
collisions can be obtained with the simple expression [30]: 



V 2 



1 



[i+k/KqY 



(3) 



where K is given by Eq. (O and Kq « 0.7 is obtained from a 
transport model calculation [31]. Thus, taking c s = 1/ v3 for 
the speed of sound of an ideal parton gas, the model has just 
two free parameters to fit to the data: the effective partonic 
cross section a gg and the elliptic flow in the hydrodynamic- 

limit \'2 lh " ■ A good agreement with Au-Au and Cu-Cu data at 

RHIC is obtained with a gg = 5.5 mb and (v^/e) = 0.22 [30] 
for a spatial density of the colliding nuclei (needed to deter- 
mine e and Ai) based on Color Glass Condensate (CGC) ini- 
tial conditions^ [32]. Namely, the elliptic flow data can be well 
reproduced with the following numerical expression 



0.22 ■ 0.45 



1 dN 1 

KoOggCs dy Ax 
dN 1 



dy A ± [mb] 



-1 



_1 dN 1 
dy Ax 

dN 1 
dy Ax[mb] 



(4) 



In this paper we take as basic assumption that the parton 
medium produced in the overlap region of p-p collisions at 
LHC energies has similar hydrodynamic properties as that in 
A-A collisions at RHIC and, thus, that the resulting anisotropic 



Hereafter, dN/dy represents the particle multiplicity per unit- 
rapidity at midrapidity, i.e. dN /dy| v= o, but for simplicity we omit the 
subindex. 

5 We note that Eq. ([4J is a "conservative" estimate of the magnitude 
of V2/e since the alternative initial Glauber matter distribution - with 
fit parameters a gg =4.3 mb and (v^/e) = 0.30 -results in a 8%-20% 
larger v'2/e ratio. 



D. d'Enterria et al: Estimates of collective hadron azimuthal anisotropy in p-p collisions at \fs = 14 TeV 



3 



flow parameter V2 follows the same eccentricity and "trans- 
verse" multiplicity scalings given by Eqs. (0 and ([3]). For sim- 
plicity we will use the same fit-parameters for a gg and ( ™ / e) 
obtained in [30], i.e. we will assume that Eq. (0]l holds too 
for p-p at the LHC. From the eccentricity £, transverse over- 
lap area Aj_ and hadron multiplicity dN/dy in p-p collisions, 
determined from a Glauber eikonal model using different phe- 
nomenological parametrisations of the proton spatial density 
distribution, we can thus estimate the expected elliptic flow pa- 
rameter, V2, as a function of the p-p "centrality" or impact pa- 
rameter b. 



at a given impact parameter b 

-b/2,y)T Pi2 (x- 
dxdyN co u^ gg (x,y;b) = 



3 Eikonal model for proton-proton collisions 

The standard procedure to determine the transverse overlap area, 
eccentricity and final multiplicities in the collision of two nu- 
clei separated by impact parameter b is based on a simple Glauber 
multi-scattering eikonal model that assumes straight-line tra- 
jectories of the colliding nucleus constituents. A recent review 
that describes the basic formalism can be found in [27]. Of- 
ten, experimentally it is more useful to introduce the "reaction 
centrality" C as a proxy for the impact-parameter b of a given 
collision, by dividing the particle production cross section into 
centrality bins Q = C\ , Ci , • • • according to some fractional in- 
terval AC of the total cross section, e.g. AC = 0.0 — 0.1 repre- 
sents the 10% most central collisions. A convenient geometri- 
cal definition of centrality is C = b 2 / (4 R 2 ), as it corresponds 
to percentiles of the total inelastic cross section in the case of 
two colliding black disks (or any other distribution with small 
diffuseness) with the same radii R, see Appendix I for details. 



Thickness and overlap functions, number of binary parton- 
parton collisions: The basic quantity of interest in a Glauber 
approach is the thickness (or profile) function of the collision 
which in the case of a nuclear reaction gives the density of 
nucleons p per unit area dxdy along the direction z separated 
from the center of the nucleus by an impact parameter b, i.e. 
TA(x,y) — A J dzp(x,y;z) for a nucleus of mass number A. By 
analogy to the nuclear case, the thickness function of a proton 
with N g partons can be written as 

T p (x,y) =N g J dzp(x,y;z), (5) 
normalised to J d 2 bT p (b) = N g . 

The overlap function of a proton-proton collision at b can be 
obtained as a convolution over the corresponding thickness func- 
tions of each proton 

T pp (b)= / dxdyT Pi i (x + b/2,y) T Pt2 (x-b/2,y) , (6) 



normalised to 



d 2 bT n 



,(b)=N 2 



For a given parton-parton cross section a gg , we can then define 
the number of binary parton-parton collisions in a p-p collision 



N C oii,gg(x,y,b) = Ggg Tp i (x- 



NcoU, gg (b) 



b/2,y) , 



*gg Tpp 



Finally, the probability density of an inelastic parton-parton in- 
teraction at impact parameter b can be defined as 



j2 pinel 

gg 

d 2 b 

j pinel 
_ gg 

db 



(b) = 



1 



Jd 2 b (l-e- 



(b) = 2%b 



or 



Jd 2 b (l-e 



(8) 



according to Poisson statistics. Note that the denominator - 
which sums over all events with at least one parton-parton in- 
teraction - is just the total inelastic p-p cross section a" 



PP 



Eccentricity and overlap area: The eccentricity of a p-p col- 
lision at impact parameter b can be obtained from the asym- 
metry ratio between the x and y "semi-axis" dimensions of the 
overlap zone, weighted by the number of parton-parton colli- 
sions at b: 



B(b) 



(y 2 -x 2 ) _ J dxdy (y 2 -x 2 ) N coU , gg {x,y;b) 



(y 2 



J dxdy (y 2 +x 2 )N co u :gg (x,y;b) 



(9) 



Other weights are possible too, see [2], but we take N co u >gg (b) 
as a natural choice - also used in more sophisticated hydrody- 
namics approaches [33] - since our picture for LHC energies is 
based on parton-parton collisions^. As a cross-check, we tested 
also a weight based on the number of participating gluons, ob- 
tained from T p (x,y) as described in [34], obtaining similar re- 
sults for z(b) for various proton densities. 

The effective transverse overlap area between the two protons 
is defined as in [30]: 



(10) 



A^b)=4nJ(x 2 )J(y 2 ) 



where the weighted averages are the same as in Eq. d9l>. We 
note that there is no commonly accepted definition of the abso- 
lute normalisation of the overlap area. Our area definition with 
maximum magnitude 4k, is four times larger than that defined 
in [23] but coincides practically with the geometrical overlap 
area of two disks with uniform two-dimensional distribution of 
density. 



Hadron multiplicity: In heavy ion collisions, the centrality de- 
pendence of the final hadron multiplicity density is found to 



6 Our eccentricity definition, Eq. ((9}, and the incomplete thermali- 
sation model one [30] are mathematically identical (modulo, a sign) in 
the case (x) : 
on r = 
and y. 



(y) = (xy) = 0, since our spatial densities depend only 
- y 2 and the N co n gg (x,y) weight is an even function of x 
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depend on a combination of the number of nucleon binary col- 
lisions and "participant" pairs [ 1 ] . Since, at LHC energies, frag- 
mentation of mini-jets produced in semi-hard parton-parton scat- 
terings is the largest contribution to midrapidity particle pro- 
duction, we will consider as a plausible operational assumption 
that the final particle multiplicity in a proton-proton collision 
follows the same impact-parameter-dependence as that of the 
number of binary parton-parton collisions N co u gg (b), i.e. that 



dN 

~dy 



(b) 



dN 

~dy 



Ncoii, gg (b) , 



(11) 



where the absolute normalisation dN /dy is chosen so as to 
reproduce the minimum-bias p-p multiplicity, namely the aver- 
age multiplicity in a proton-proton collision integrated over all 
impact-parameters 



dN m 
dy 



dN 
dy 



d 2 bN co ii, gg (b) 



2 pinel 

gg 



d z P: 



d 2 b 



(b), 



(12) 



which, at midrapidity at LHC energies, is expected to be 
dN MB /dy fa 10 according to the PYTHIA [10] or PHOJET [35] 
MCs based on various extrapolations from lower energy data. 



3.1 Proton spatial transverse densities 

The matter distribution in transverse space of the colliding pro- 
tons, p (x, y) determines the probability to have a given num- 
ber of parton-parton interactions at impact parameter b. In or- 
der to quantify the proton profile function, Eq. ©, one often 
assumes a spherically symmetric distribution of matter in the 
Breit system, p{x)d i x — p(r)d^x, and, for simplicity, one takes 
the same spatial distribution for all parton species in the proton 
(valence and sea quarks, gluons) and momenta, and neglects 
any correlations among them. The spatial proton density has 
been parameterised with different distributions in the literature 
(see e.g. [8, 10]) such as: 

- Hard sphere: The simplest model is to consider that the pro- 
ton has a spherical form with uniform parton density of ra- 
dius R: 

pW = 47^ e < r -*> (13) 

The associated root-mean-square (rms) charge radius is 

3/5/?. Of course this is an unrealis- 



Rrms — V V ) — 

tic approximation of the proton shapfl but we keep it as an 
"extreme" case since the overlap area of two such distribu- 
tions has very large eccentricities. 
- Exponential: Matter can be distributed in the proton accord- 
ing to its charge form factor, which is well represented by 
an exponential expression in coordinate space: 



Pto = 



1 



,-r/R 



(14) 



reproducing to a large extent the spatial distribution of the 
valence quarks, with rms radius R rms = y/T2R. 



1 E.g. in particular at higher virtualities, such a profile violently dis- 
agrees with the small-/ (where t is the squared momentum exchanged) 
J/\if data at HERA [36]. 



Fermi distribution: The standard spatial density for nuclei 
is the Fermi-Dirac (or Woods-Saxon) distribution with ra- 
dius R and surface thickness a 



P(r) 



Po 



,{r-R)/a + 1 ' 



(15) 



where po is a normalisation constant so that J d 3 rp(r) = 1. 
The associated rms radius is R rms = y/3/5R in the limit 
a/R — > and R rms = \f\2a for a/R — > °°. For a/R = 0.2, 
R rm w 1.07/?. 

Gaussian: A simple Gaussian ansatz, although not very re- 
alistic, makes the subsequent calculations especially trans- 
parent and, therefore, is often used in the literature: 



1 



( V2kR) : 



exp 



2R 2 j 



(16) 



The corresponding rms radius is given by R rms = \/2>R. 
- Double-Gaussian: This corresponds to a distribution with a 
small core region of radius «2 containing a fraction p of the 
total hadronic matter, embedded in a larger region of radius 
a\: 



i-P 



exp < ? 



P 

— exp 

a 2 



(17) 



Although the Gaussian and double-Gaussian are popular pro- 
ton density choices (e.g. in the PYTHIA MC), we will not con- 
sider them hereafter because mathematically the convolution 
of two such distributions results respectively in an exactly null 
or very small eccentricity, i.e. they cannot generate any elliptic 
flow according to our ansatz, Eq. (0). 

Table Q] lists the input parameters used for three of the five 
proton densities mentioned above. The chosen radius R of the 
hard-sphere parametrisation is consistent with electron-proton 
scattering fits [37] and with diffractive results at HERA (more 
sensitive to the gluon density) [38]. The parameters of the ex- 
ponential (R) and Fermi-I (R and a) distributions are those ob- 
tained by the CDF collaboration in their data analysis of dou- 
ble parton scatterings in p-p collisions at y/s =1.8 TeV [8]. In 
addition, we quote the parameters of the proton Fermi-II dis- 
tribution obtained from a Fourier transformation of the energy- 
dependent p-p elastic amplitude as discussed in [39]. Such an 
approach aims at taking into account the effective growth of 
the proton size due to the larger transverse spread ("diffusion") 
of partons in the proton for increasing collision energy. Within 
this framework, the radius of the strong interactions is expected 
to be a factor 1.5 larger at the LHC compared to fixed-target 
energies. From unitarity (i.e. the optical theorem), the squared 
cm. energy (s) dependence of the inelastic (total minus elas- 
tic) p-p cross section can be written as an integral over impact 
parameters [40]: 

=<W"<W = / d 2 b [2Rer pp (s;b)-\r pp (s;b)\% 



(18) 

where F pp (s;b) is the profile function of the elastic amplitude 
and the integrand represents the distribution of the cross section 
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Table 1. Top rows: Radius R and diffusivity a parameters for various proton spatial densities, Eqs. J X 3 b - l |15l >, and effective number of gluons 
N g in the proton, considered in this work. The three bottom rows show the maximum eccentricity and overlap area and inelastic cross section 
derived for p-p collisions at y 7 ? = 14 TeV from each set of parameters within our Glauber approach. The last column shows the corresponding 
values for typical Au-Au collisions at RHIC energies. 





Hard sphere 


Exponential 


Fermi-I 


Fermi-II 


Fermi Au-Au 




refs. [37,38] 


ref. [8] 


ref. [8] 


ref. [39] 


ref. [28] 


radius R (fm) 


0.89 


0.20 


0.56 


1.05 


6.36 


diffusivity a (fm) 






0.112 


0.29 


0.54 


effective number of partons N g 


17 


9 


17 


4 


197 (nucleons) 


rms radius R rms (fm) 


0.69 


0.70 


0.60 


1.34 


5.32 


diffusivity/radius a/R 






0.2 


0.28 


0.085 


effective partons/fm 3 N g p(r = fm) 


5.8 


44.8 


16.7 


0.48 


0.17 (nucleons/fm 3 ) 


p-p eccentricity e max at b ~ 2R rms 


1. 


-0.4 


0.07 


-0.05 


0.4 


p-p overlap areaA^ at b = fm (fm 2 ) 


1.6 


0.58 


0.85 


3.8 


85. 


p-p inelastic cross section o™p' (mb) 


80. 


79.6 


78.5 


90. 


7110 



for generic inelastic collisions (i.e., summed over all inelastic 
final states) over impact parameters. The ^-dependence of the 
inelastic probability defined as 



j1 pinel 

88 

d 2 b 



(s;b) 



2Rer pp (s;b)-\r pp (s;b)\ 



vine 1 1 
J PP 



(19) 



can be fitted to the available p-p elastic data at high-energies [41] 
and extrapolated [39] to LHC energiefl By fitting P'™ l (s;b) to 
Eq. dl5t . we obtain the R and a parameters for the Fermi-II spa- 
tial density quoted in Table Q] 

For each one of the density parametrisations considered, we 
obtain a proton root-mean square (rms) charge radius not far 
from that obtained from the world-data on e-p elastic scattering 
at low momentum transfer (Q < 4 fm~ '): R rms w 0.89 fm [42], 
except for the Fermi-II case which, as aforementioned, yields 
an effective proton radius twice larger at the LHC than ob- 
served at low energies. For comparison purposes with p-p, we 
add also in the last column of Table Q] the results for Au-Au 
collisions at ^fs^, = 200 GeV with the standard Au Fermi dis- 
tribution parameters [28]. 

The parton-number normalisation N g of the proton thick- 
ness function, Eq. (0, can be obtained by requiring that the 
proton-proton inelastic cross section obtained from the p-p over- 
lap function, Eq. (0, in the eikonal approximation 



-inel 

PP 



d 2 b 



(20) 



We note that the parton-parton cross section, G gg = 7.8 mb, ob- 
tained from the fit of [39] is not far from the a gg = 5.5 mb obtained 
within the incomplete thermalisation model, Eq. ©. Both should be 
interpreted as effective semihard partonic cross-sections, larger than 
the standard perturbative gluon-gluon cross-section which at LO is 
a gg = 9/2na 2 /p 2 fs 2 mb for oc, = 0.5 at a pr-cutoff of order 
H = 1 GeV. 



is in the range expected at LHC energies, a 1 "* fts 80 mb [10, 
35], for our considered effective parton-parton cross section of 
G gg = 5.5 mb. The corresponding number of partons normali- 
sation N g is shown in the Table [T]for each proton matter distri- 
bution. We also tabulate the product N g p(r = fm) indicating 
the parton density in the centre of the proton. For comparison, 
we also show the density in the center of the Au nucleus (with 
A = 197 nucleons), which is equal to the well-known value 
Ap&(r = fm) =0.17 nucleon/fm 3 . The corresponding inelas- 
tic Au-Au cross section at ^fs^, = 200 GeV, GauAu = 7110 mb, 
obtained from Eq. < f20b with a NN = 42 mb, is larger as expected 
than that of two black absorptive hard disks, GauAu = 5080 mb, 
as the Fermi profile for both nuclei has long tails. 



3.2 Number of binary collisions and parton-parton 
inelastic probability 

In any eikonal approach, once the spatial densities of the two 
colliding objects are known, one can derive any geometrical 
quantity related to their interaction from the overlap function, 
Eq. (|7l), or equivalently from the number of binary collisions 
of their constituent particles, Eq. (0. In the p-p case, the num- 
ber of parton-parton collisions N co u,gg as a function of scaled 
impact parameteiQ b/2R and centrality C = b 2 / (4 R 2 ) for the 
different proton spatial densities considered in this work are 
shown in the left and right plots of Figure Q] respectively. The 
plot as a function of centrality is shown as an histogram with 
binning AC = 0.2 as could be obtained experimentally by di- 
viding the data in 20% fractions of the measured cross section 
(e.g. based on an observable monotonically increasing with im- 
pact parameter, such as the event particle multiplicity). 

The exponential and Fermi distributions feature a long tail 
well beyond twice the radii R, whereas the hard-sphere parametri- 
sation dives off quickly when b approaches 2R. For the most 



9 For the exponential density, for which R and R rms are quite differ- 
ent, we use b/(2R rms ). 
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b/2R C 

Fig. 1. Number of binary parton-parton collisions in p-p collisions at ^/s = 14 TeV as a function of scaled impact parameter b/2R (left) and 
centrality C (right) for the different proton density distributions considered in this work (Table Q}. For comparison, the results for Au-Au at 
RHIC energies are shown as a dotted line. 




Fig. 2. Probability density (normalised to unity) of inelastic parton scatterings in p-p collisions at ^/s = 14 TeV as a function of scaled impact 
parameter b/2R (left) and centrality C (right) for the different proton density distributions considered in this work (Table[TJ. For comparison, 
the results for Au-Au at RHIC energies are shown as a dotted line. 



central p-p collisions (b = fm), the hard-sphere and Fermi - 
I distributions yield about 100 parton-parton interactions, the 
exponential gives about 30 collisions, and the Fermi-II distri- 
bution, which has a larger proton size and thus a more "dilute" 
parton density, just 2 collisions. For comparison, the number of 
nucleon-nucleon collisions in Au-Au at RHIC shows a steeper 
centrality-dependence than the Fermi densities for p-p colli- 
sions because the nuclear distribution has a relatively smaller 
edge diffuseness (a/R = 0.085) than the proton one. 



The other key quantity of our approach is the impact-parame- 
ter dependence of the midrapidity particle multiplicity produced 
in p-p collisions at = 14 TeV. Within our framework, we 
estimate dN/dy(b) via Eqs. (fTTT) and ST2\ . The latter depends 
on the parton-parton inelastic collision probability dP' g " g el /db, 
Eq. (8), which is shown as a function of scaled impact-parameter 
(left) and centrality (right) in Fig.|2]for the various proton den- 
sities considered in this work. All probabilities are normalised 
to unity to facilitate the comparison. 
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4 Eikonal model results for p-p collisions at the 
LHC 

Based on the eikonal model formalism and proton density parame 
trisations discussed in the previous section, we derive here all 
basic quantities needed later to estimate the expected amount of 
elliptic flow in proton-proton collisions at y/s = 14 TeV: eccen- 
tricity, transverse overlap area and particle multiplicity. Again, 
in order to compare the results for the various proton densi- 
ties considered in Table Q] (with varying radii R), we present 
our results as a function of the dimensionless scaled impact- 
parameter bj (2R) and centrality C = b 2 / (AR 2 ), see Appendix! 

Eccentricity: Figure|3]shows the computed eccentricity in p-p 
at y/s = 14 TeV as a function of the scaled impact-parameter 
(left) and centrality (right) of the collision. The eccentricity is 
increasingly large and positive for the hard-sphere parton den- 
sity (up to e = 1 for b = 2R) and increasingly negative (down 
to e = —0.4 for b « 2R) for the exponential spatial density. 
For the Fermi distributions with diffuse edges, the eccentric- 
ity is small (e < 0.1) and changes its sign for very large im- 
pact parameters. The larger the ratio of diffuseness-to-radius 
(a/R = 0.28, 0.2, and 0.0085 for the Fermi-II, Fermi-I, and 
Fermi-AM-Ai< respectively) is, the smaller the impact-parameter 
at which e changes sign. In any case, the small magnitude of 
the p-p eccentricities for the Fermi-I (about five times smaller 
than Fermi-AuAu) and Fermi-II (very close to zero) profiles 
and the large impact-parameters where the sign change takes 
place, likely preclude to see any such effect in the elliptic flow 
data. Also, as aforementioned, it is easy to show that the over- 
lap of two Gaussian distributions for the colliding protons, has 
zero eccentricity at any impact parameter. 

Overlap area: Figure|4]shows the p-p transverse overlap area 
as a function of scaled impact parameter. The left plot shows 
the absolute A± (b) value, and the right plot presents the overlap- 
area normalised to the maximum value for a "head-on" colli- 
sion at zero impact-parameter, A± (b)/A± (b — fm) . The over- 
lap area for central p-p collisions is largest for the Fermi-II 
density (Aj_ rj 4 fm 2 ) and smallest for the exponential distri- 
bution (Ax ~ 0.6 fm 2 ). Again, as for the eccentricity, the hard- 
sphere and exponential distribution have opposite behaviours 
as a function of centrality: A±(b) decreases (increases) for the 
former (for the latter). On the other hand, there is only a small 
dependence of the transverse area on impact-parameter for the 
overlap of two Fermi densities with diffuse edges. Compared 
to Au-Au, the obtained values of A± in p-p collisions, for the 
Fermi and hard-sphere densities (Aj_ w 1 fm 2 ), are about two 
orders of magnitude smaller. 

Hadron multiplicity: The third ingredient of the incomplete 
thermalisation model, Eq. (01, is the particle multiplicity per 
unit rapidity at midrapidity, dN/dy, obtained via Eqs. ( fTTT i - 
(TT2l) . The results are shown in Fig. [5] for a valu43 of dNo/dy 

10 Note that, given our absolute normalisation based on dN MB /dy, the 
actual N g and a gg values do not play practically any role in the final 
determination of the hadron multiplicity. 



chosen in Eq. (fT2l) so as to reproduce the expected dN MB /dy — 
10 multiplicity in minimum-bias p-p collisions at the LHC. All 
proton densities yield similar values of dN/dy in a wide range 
■of centralities C in p-p collisions. The 20% most central colli- 
sions result in 2-3 times larger multiplicities than in minimum 
bias p-p (or Au-Au) collisions for all considered densities. For 
the comparison Au-Au data, we take dN/dy = 600. 

As we see from our elliptic-flow "ansatz", Eq. @, the V2 
parameter depends on the particle multiplicity but only after 
normalisation by the size of the system. Figure [6] shows the 
transverse density multiplicity, namely the ratio of the hadron 
multiplicity over the transverse overlap area (dN/dy) /A±, for 
different forms of the proton matter distribution. We find ratios 
(dN /dy) /A± w 2.5 - 4.5 mb _1 for all distributions in cen- 
tral p-p collisions similar to those found in Au-Au collisions 
at RHIC except for the dilute Fermi-II density, which features 
a much lower (dN /dy)/Aj_ w 0.5 mb~' value. Interestingly, 
although the particle multiplicity is about two orders of magni- 
tude larger in nuclear than in proton collisions, the overlap area 
is roughly smaller by the same amount (see Fig. |4j and, there- 
fore, the multiplicity density is not very dissimilar in p-p and 
A-A collisions: (dN/dy)/Aj_ wO.l -4.5 mb _1 , for all densities 
within a large interval of centralities. 

5 Elliptic flow estimates in p-p collisions at the 
LHC 

From the values of eccentricity e (Fig. [3), overlap area A± 
(Fig.|4]i and hadron multiplicity dN/dy (Fig. [5), as a function of 
impact-parameter b determined in the previous chapter, and us- 
ing the Eq. dU of the incomplete thermalisation model, we can 
finally obtain the value of the integrated elliptic flow parameter 
V2 for each one of the four different proton densities. By con- 
struction, \'2{b) will follow the general behaviour of e(b). Fig- 
ure|7]shows the integrated \'i as a function of transverse particle 
density expected at midrapidity in p-p collisions at LHC ener- 
gies. In Fig. [8] we show vj versus centrality (left) and versus 
particle multiplicity (right) normalised by the maximum value 
in the most central events, N(C) /N max , which is a proxy for the 
reaction centrality C often used experimentally. 

Several features are noticeable in both figures. First, the 
hard-sphere distribution features a maximum value of V2 « 0. 1 , 
almost twice higher than in Au-Au collisions and much larger 
than the rest of p-p elliptic flow parameters. The fact that the 
overlap of two hard-spheres with infinitely sharp edges (see 
Fig. [3]l yields artificially large eccentricities is a well-known 
fact [2]. We plot this result just to delimit the theoretically 
maximum amount of V2 that could be generated within our 
approach. Second, the maximum v>2 obtained with the Fermi- 
I profile is about 1.5% and vanishingly small for the dilute 
Fermi-II distribution. Third, exponential transverse densities 
for the proton, feature negative p-p eccentricities with mini- 
mum elliptic flow parameters of the order of —3.5%. A nega- 
tive integrated V2 would indicate that the anisotropy of hadron 
emission is not in-plane (as for a positive v'2) but out-of-plane, 
i.e. perpendicular to the reaction plane. 
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Fig. 3. Eccentricity e in p-p collisions at y/s = 14 TeV as a function of scaled impact parameter b/2R (left) and centrality C (right) for the 
different proton density distributions considered in this work (TableQJ. For comparison, the results for Au-Au at RHIC energies are shown as a 
dotted line. 




b/2R 



Fig. 4. Effective overlap area Aj_ in p-p collisions at y^j = 14 TeV as a function of scaled impact parameter b/2R for the different proton 
density distributions considered in this work (TableQJ. For comparison, the results for Au-Au at RHIC energies are shown as a dotted line. The 
left plot shows the absolute value of A ± (for clarity, the Fermi-II and Au-Au curves are scaled by factors of 1/4 and 1/100 respectively). The 
right plot shows the area normalised to the value for central collisions, A± (b) /Ax (b = fm) . 



Setting aside the obtained with the unrealistic hard-sphere 
profile, the maximum absolute values of the integrated elliptic 
flow in p-p collisions in the range l%-3% are consistent with 
those obtained in a recent hydrodynamics study [17] - which 
uses a spatial density based on a parametrisation of the Fourier 
transform of the proton electromagnetic form-factor [43] - which 
predicts also v>2 < 0.035 (the maximum value obtained when 
the produced strongly interacting matter has vanishing viscos- 
ity) at 14 TeV. Also percolation models [15] predict integrated 



V2 of the order of 2%-3%. Although small, such v'2 values are 
comparable to those measured in nucleus-nucleus reactions at 
cm. energies i/s^ « 5 GeV at the BNL Alternating Gradient 
Synchrotron (AGS) [44] or y/s^ « 17 GeV at the CERN Super 
Proton Synchrotron (SPS) [26] and thus could be in principle 
measurable. 

The experimental measurement of j V2 1 signals of maximum 
magnitude £?(3%), in p-p collisions at the LHC will be cer- 
tainly challenging. On the one hand, one will have to deal with 
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Fig. 5. Absolute (left) and normalised (right) particle multiplicity at midrapidity as a function of centrality C in p-p collisions at = 14 TeV 
for the different proton density distributions considered in this work (Table [TJ. For comparison, the results for Au-Au at RHIC energies are 
shown as a dotted line. 
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Fig. 6. Particle multiplicity density at y = per unit of transverse overlap area Ax in p-p collisions at yfs = 14 TeV as a function of centrality 
C for the different proton density distributions considered in this work (TableQJ. For comparison, the results for Au-Au at RHIC energies are 
shown as a dotted line. 



all standard "non-flow" effects - i.e. azimuthal correlations not 
associated with parton collective rescattering - generated by 
jets, resonance decays, Hanbury-Brown Twiss (HBT) short- 
range correlations, or (depending on the technique used to de- 
termine v'2) momentum conservation [2]. In particular, the large 
minijet production cross section at the LHC generates already 
strong (back-to-back) azimuthal anisotropics [45,46]. In ad- 
dition, part of the azimuthal correlations can come, not from 
final-state rescatterings but from parton correlations in the ini- 
tial state [16,47]. On the other hand, the relatively small magni- 
tude of the p-p eccentricity and overlap area will result in large 



fluctuations of the v'2 parameter even for a same given central- 
ity. Likewise, event-by-event fluctuations on the produced par- 
ticle multiplicities will also complicate their correlation to the 
measured v'2 values, as commonly presented in plots like that 
in Fig. [8] (right). In any case, there exist detailed methods of el- 
liptic flow analyses developed for nuclear collisions [22,23,25, 
48-51] that can be tested and adapted to p-p collisions at the 
LHC. We leave for a future paper the detailed study of the ex- 
perimental feasibility of the measurement. As aforementioned, 
previous studies in heavy ion collisions at lower energies have 
clearly demonstrated the possibility to measure V2 signals with 
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Fig. 7. Integrated elliptic flow parameter as function of transverse particle multiplicity, (dN /dy)/A^, expected at midrapidity in p-p 
collisions at *Js = 14 TeV for the different proton density distributions considered in this work (Table[T}. For comparison, the Au-Au results are 
shown as a dotted line. 




Fig. 8. Integrated elliptic flow V2 parameter as function of centrality (left panel) and of normalised particle multiplicity (right panel) at 
midrapidity in p-p collisions at ^/s = 14 TeV for the different proton density distributions considered in this work (Table[TJi. For comparison, 
the V2 for Au-Au at RHIC energies is shown as a dotted line. 



absolute magnitudes of order l%-3% and the ALICE [52], AT- 
LAS [53] and CMS [54] experiments have proven detector ca- 
pabilities to carry out such measurements. 



6 Summary and conclusions 

We have studied the possibility to observe a collective expan- 
sion signal - in the form of an azimuthal anisotropy of particle 
production with respect to the reaction plane - due to multi- 
parton interactions in proton-proton collisions at the LHC. Our 
working assumption has been that any possible azimuthal aniso- 



tropy due to collective flow in p-p collisions should follow the 
same eccentricity, overlap-area and particle multiplicity depen- 
dences observed in the strongly interacting matter formed in 
high-energy heavy ion collisions. Using a simple eikonal model 
for multiparton scatterings, we have tested various proton den- 
sity distributions proposed in the literature and obtained the 
corresponding eccentricities, transverse areas and hadron mul- 
tiplicities as a function of the impact parameter for p-p col- 
lisions at y/s = 14 TeV. The transverse overlap area A± and 
the final particle multiplicity at midrapidity dN/dy are about 
two orders of magnitude smaller in p-p compared to Au-Au 
collisions. Since the elliptic flow roughly depends on the nor- 
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malised (dN/dy)/A± ratio, we thus expect any elliptic flow in 
the p-p systems to be mostly driven by the eccentricity (if any) 
of the system produced in the collision. 

Our first finding is that, popular proton matter profiles such 
as Gaussian or double-Gaussian result in vanishing p-p eccen- 
tricities and, thus, cannot generate any final-state elliptic flow 
within our approach. Unphysical sharp edge (hard-sphere) pro- 
files result in integrated elliptic flow parameters, v>2 ~ 10%, 
almost twice larger than found in A-A at RHIC. More realis- 
tic Fermi-Dirac distributions with a diffuse proton edge, yield 
maximum V2 values of the order of 1.5%. If the Fermi density 
is too dilute, as in our considered Fermi-II case that takes into 
account an effective growth of the proton size due to the trans- 
verse spread of its partons at high energies, the generated V2 
will be virtually null. Lastly, exponential proton profiles that 
reproduce the proton charge form-factor (i.e. the spatial distri- 
bution of its valence quarks) would result in negative, i.e. out- 
of-plane, integrated elliptic flows with minimum values of the 
order of —3%. 

All in all, our work demonstrates that the study of hadron 
anisotropies with respect to the reaction plane in p-p collisions 
at LHC energies, can provide important information on the pro- 
ton shape and structure at moderate virtualities 0(0.7 GeV). 
Although previous analyses with heavy ions at much lower en- 
ergies have indeed measured integrated V2 of a few percents, the 
experimental extraction of such a signal in p-p collisions will 
be challenging given the expected large non-flow azimuthal 
correlations that can mask the signal. Despite these difficul- 
ties, the absence or presence of elliptic flow in the data and its 
dependence on the "centrality" of the collision, will nonethe- 
less put strong constraints on the density profile of the proton 
at very high energies. 
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Appendix I. 

We recall here the basic transformations from impact parameter 
to centrality for a general two-dimensional probability density 
F(b) which depends on the impact parameter vector b, 



d 2 P 
d 2 !? 



d 2 P 1 dP 

= =F(b), 

bdbdv? 2%bdb v ; ' 

d 2 P 



f 2 d p f , 
with normalisation J d b—^g — J db(2nbF{bj) = 1 



The corresponding normalised one-dimensional density for the 
dimensionless impact parameter b' —b/ (2R), which we have 
used in this work to easily compare the results obtained with 
proton densities with different radii R, is 

dP dP 

— = 2%bF(b) -> — =A%RbF(b), 

db db' 

For our alternative centrality parameter defined as C — b 2 / (4-R 2 ), 
i.e. db = 2R 2 /b dC, we have 

dP 2R 2 dP 2R 2 . . 7 , s 

—2%bF(b)=4nR 2 F(b) , 



normalised to / dC 



dC b db 
dP 



dC 



b 2R 2 dP 

— ~db 

2R 2 b db 



db- 



dP 

) 

db 



1. 



Such a definition of C is convenient as it corresponds to a given 
fraction of the inelastic cross section. Indeed, for two colliding 
black disks with radii R\ and R2 the fraction of the inelastic 
cross section is 

Ab 2 

b < (Ri +R 2 ) , 



Ao,„ e / 



2nbAb 



Oinel n{R l +R 2 ) 2 (Rl+Rl) 2 

and thus for two identical disks R = R\= R2, we have 



Agingj _ A / b 2 



AC. 



The exact definition of centrality c as a fraction of the inelastic 
cross section is equal to [55] 



'geom\ u cJ 



2%lbbdb{V - g-°« r w>(*)) 



•zinel 



In Fig. [9] we compare the exact centrality c to the centrality 
C used in the present work. Let us consider, as an example of 
b-io-C transformation, the case of the probability of inelastic 
parton scatterings, presented in Section [3] Such a distribution, 
which as a function of impact parameter and centrality reads: 



//if/ 



dP" 

88 
db 



- 2%bF(b) 



2nb(l - e~ a «* r wW) 

J d 2 b(l - e~ a ss Tp P (b)^ 



and 



—$$-=4nR 2 F(b)= " { — /, b = 2RVc 

dC y ' J d 2 b(l - e- G ss T PP (b)y 

is used in our calculations to determine the multiplicity density 



d 3 N 
dyd 2 b 



(b) 



dN M , 
dy 



Ncoii, sg {b) -3%-(b) 

a2 pinel 

fd 2 bN coll , gg (b)-^(b) 



The multiplicity for different centrality bins is: 

dy Jc 



Q+AC ^jy 

c k dydC 



The dimensionless multiplicity and probability in a given impact- 
parameter and centrality bin are 

dN dN 
N(b) = -Ab, or N(C) = -AC, 



(b) 



j pinel 

db 



or i*f(C) 



dP'"f' 
—^AC. 
dC 



12 



D. d'Enterria et ah; Estimates of collective hadron azimuthal anisotropy in p-p collisions at \fs = 14 TeV 



E 



E 1 

o 

to 

0*0.9 
~~ 0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



pp Hard-sphere 
pp Fermi-I 
pp Fermi-ll 
pp Exponential 
AuAu Fermi 



" I 



1 



_J 



J— ■ 




0.2 



0.4 



0.6 



0.8 



1.2 



1.4 

C 



Fig. 9. Comparison of the exact centrality c, Eq. J2 1 b , to the centrality C = b 2 /(4R 2 ) used in the present work. 
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